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Abstract 

The Eulerian triangle is a classical array of combinatorial numbers 
defined by a linear recursion. The associated boundary problem asks 
one to find all extreme nonnegative solutions to a dual recursion. Ex- 
ploiting connections with random permutations and Markov chains we 
show that the boundary is discrete and explicitly identify its elements. 



1 Introduction and main results 

The Eulerian triangle (see, e.g., [HI Section 6.2]) is the infinite array of Eu- 
lerian numbers 1 

n \ 

k ) (0< k<n-l, n= 1,2,...), 
considered together with the defining them recursion 

;h* +i ><Y) + <»-*>(;:|) « 

and the boundary conditions 

\=1, and (^) = for n<0 or fc>n-l. (2) 
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The first six rows of the triangle (n = 1, . . . , 6) are 
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We are interested in nonnegative solutions 

V = (V nk : < k < n - 1 , n = 1, 2, . . . ) 

to the dual or backward recursion 

Vnk = {k + l)K+i,fc + (n - A;)K + i )fc+ i, (3) 

subject to the only normalization condition V\q = 1. In contrast to the 
dual recursion has multiple solutions which comprise a convex set V. We 
denote by ex(V) the set of extreme solutions and call it the boundary. 

One principal result of this paper gives a natural parametrization of the 
boundary: 

Theorem 1. Every extreme solution W = (W n k) £ ex(V) is uniquely de- 
termined by the parameter 9 := W20, which assumes values in the following 
subset of the unit interval 

\2 x+1' J \2J \2 x+1' J 

The correspondence between ex(V) and is a homeomorphism. 

Thus, the parameter set 9 is composed of two sequences and their sole ac- 
cumulation point ~. An obvious symmetry of 9 about | corresponds to the 
symmetry of (HJ) under the substitution (n, k) — > (n, n — 1 — k). 

Theorem^distinguishes the Eulerian triangle from other classical number 
triangles, whose boundaries can be also identified with subsets of the unit in- 
terval. It is well known that the boundary of the Pascal triangle is [0, 1]; a re- 
sult equivalent to both de Finetti's theorem on exchangeable trials and Haus- 
dorff's characterization of moment sequences [TJ |3J H21 122] ■ The boundary 



2 



of the g-Pascal triangle of g-binomial coefficients is {0} U {. . . , q~ 3 , q~ 2 , g -1 } 
for q > 1, see IT7] . For a parametric family of the generalized Stirling 
triangles considered in [7j the boundary was shown discrete (as for the q- 
Pascal triangle) for some values of the parameter and coinciding with [0, 1] 
for other. For instance, the boundary is discrete for the triangle of Stirling 
numbers of the second kind, and it is continuous for the triangle of signless 
Stirling numbers of the first kind. Some sporadic earlier results on a wider 
family of generalized Stirling triangles are found in [TU] . 

The extreme solution corresponding to a particular value of the parameter 
9 e 6 will be denoted W{9) = (W nk (9)). Our second principal result gives 
explicit formulas for these solutions: 

Theorem 2. For 9 = 2( ^| f 2 1) > \ with x = 0, 1, 2, . . . we have 

W nk (9) = ( H + H ~ k ) /(H+l) n (n = l,2,...;0<A;<n-l), (4) 



n 

whereas for 9 = 2 ^ +l ^ < \ with x = 0, 1, 2, . . . we have 

W nk {9)= ( X + ^ + 1 ) /Wir (n = l l 2,...;0<*<n-l). (5) 



Finally, the solution corresponding to 9 = \ 



is 



W nk (\) = ^ (n=l,2,...; 0< k<n-l). (6) 



It is seen that for 9 > |, W nk {9) is for k > x, and for 9 < | it is for 
k < n — 1 — x. Notice also the symmetry 

W nk {9) = W n>rM {l-9). 

All three formulas of Theorem |21 can also be written in a unified way 

— k+n— 1 

Wn*(0) = ^ J] ( 1 + ^)> ( 7 ) 

i=— k 

where 9' := 29 — 1 ranges over the set 

0' = (-1 —- —- —- i i i H 
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Theorems n and 121 are our main results. They imply a simple description 
of the whole set of nonnegative normalized solutions to (jSJ) • 

Corollary 3. Each solution V G V can be uniquely represented as a convex 
combination 

with p a probability distribution on 0. 

The rest of the paper is organized as follows. In Section |2] we encode the 
Eulerian triangle into a graded graph £ , which allows to embed our boundary 
problem in the general framework developed in PUJ EU E3 E] • In particular, 
we recall a construction of the Martin boundary, which is our major techni- 
cal tool. In Section |3] we explain connections to random permutations and 
random total orders onN = {1,2,...} with a sufficiency property. In Section 
EJthe extreme solutions are related to random sorting algorithms. The main 
results are proved in Sections El and El Some connections and extensions are 
discussed in the last section. 

Boundary problems of combinatorial nature arise within a variety of 
mathematical disciplines like probability, numerical analysis, algebra and rep- 
resentation theory. Our own view on the subject has two major sources: the 
asymptotic representation theory of symmetric groups initiated by Kerov and 
Vershik '24\ and Kingman's theory of exchangeable partitions summarized 
in lectures by Pitman [22 J. 

2 The graph S and the Martin boundary 

The graph £ is an infinite directed graded graph with the vertex set 

{(n,k); n = 1, 2, . . . , < k < n - 1} 

and multiple edges. The nth level of £ , denoted £ n , is the set of n vertices 
(n, 0), . . . , (n, n — 1). The level £\ has a single vertex (1, 0), which is the root 
of £. For generic vertex (n, k) the outgoing edges link (n, k) to two vertices 
(n + 1, k) and (n + 1, k + 1) by k + 1 and n — k directed edges, respectively. 
We say that two vertices connected by an edge are adjacent with each other. 
A path in £ (finite or infinite) starting at (n, k n ) is a sequence of edges 
linking adjacent vertices (n,k n ),(n + l,k n+ i), . . . on consecutive levels. A 
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path starting at the root (1,0) is called standard. The edge multiplicities in 
£ are selected to match with the coefficients in the recursion from which 
it is clear that the Eulerian number (") is the dimension 2 of the vertex 
(n, k), meaning the number of standard paths in £ which terminate at (n, k). 

The interpretation of dual recursion Q requires some concepts of proba- 
bility theory. Let us consider £ as a state space of some Markov process whose 
time parameter n runs in the reverse direction >n—>n — 1— 1, 
a possible state at time n is a vertex in £ n , and the transition probabilities 
are given by 

Prob{(n, fc) - (n - 1, k)} = (k + 1) (8) 

\fc/ 

/n-l\ 

Prob{(n, k) - (n - 1, k - 1)} = (n - k)-j^- . (9) 

\ fe/ 

The basic relations (JTJ) and ((21) translate as the rule of total probability 



Prob{(n, k) -> (n - 1, fc)} + Prob{(n, jfe) -> (n - 1, A; - 1)} = 1, (10) 

and imply that at consecutive times the process must reside in adjacent 
vertices of £. 



Now let V = (V nk ) e V. Setting 



V nk = Q)v nk , (11) 



the recursion Q can be rewritten as 

V nk = Prob{(n + 1, k) -> (n, fc)}y n+ i )fc 

+ Prob{(n + 1, fc + 1) -> (n, fc)}K+i, fc+ i 

Lemma 4. VKe /lave 

n-l 



(12) 



53^ = 1, n=l,2,.... (13) 



fc=0 



2 The dimensions alone do not determine £ unambiguously, since the Eulerian numbers 
satisfy many other recursions different from Q . By our definition of the Eulerian triangle 
we were careful to say that the array was considered together with the recursion Q . 
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Proof. Indeed, from (JTUJ) and (fT2j) . the quantity Ylk=o does not depend 
on n. Since it equals 1 for n = 1 due to the normalization condition, the 
same holds for all n. □ 

Thus, the vector (V n o, . . . , V n>n -\) is a probability distribution on £ n for each 
n, and this family of distributions is consistent with the transition probabil- 
ities (JH1) and 0- It follows that V determines the law of a Markov chain by 
the virtue of (fTTj) . The boundary problem acquires therefore the following 
meaning: 

• describe all possible probability laws for a Markov chain on S, whose 
transition probabilities are given by (jHJ) and 0. 

If we required © to only hold for n restricted to some finite range 1 < 
n < N, the analogous boundary problem were rather simple. For, each 
truncated solution 

(V nk ; 0< A; < ri-l, n = l,...,N) 

is uniquely determined by the last row (Vjyo, . . . , Vn,n-i)- Equivalently, the 
corresponding Markov chain with time parameter n ranging from N to 1 
is determined by the initial distribution (Vno, ■ ■ ■ ,Vn,n-i), which can be 
selected arbitrarily within the set of all probability distributions on Sn- The 
set of (nonnegative, normalized) solutions to such a truncated recursion 
is therefore the convex hull of the arrays V N>C (0 < x < N — 1) which have 
the iVth row 

Vg?=8 kh l ( N ) (0<k<N-l), 



where 6^ is the Kronecker symbol. The set is a (N — l)-dimensional 
simplex with extreme elements 

ex(V (A °) = {V Nit : x = 0,...,iV-l}. 

The probability law V N>C corresponding to V N * rules a Markov chain which 
starts in state (N, x) G £n at time N, hence the boundary of the iV- 
truncated triangle can be identified with £ N . 

For the infinite recursion the problem is much more complicated because 
there is no obvious analogue of the "last row" which would provide an initial 



6 



condition for (J3J). A common recipe to obtain all solutions is the following. 
Extend each V N * to a function on the whole set of vertices of £ by setting 
V^fc C = for n > N. Define the Martin boundary 3 of £ as 



where the bar means the closure in the topology of pointwise convergence of 
functions on the set of vertices. Plainly, dV is the set of solutions which may 
be obtained from truncated solutions V Nit by fixing some limiting regime 
for x = *c(N), as iV — > oo, to secure convergence of V^T f° r eacn (n,k). 
Obviously, each such limit is indeed a solution to ©, hence dV is a subset 



By some well known general theory (see Ch. 1, §1] and also PQEUHl 
IT2] ) the Martin boundary dV contains the boundary ex(V) (for this reason 
dV is sometimes called the maximal boundary). In our proof of Theorems^ 
and 121 we shall determine the Martin boundary dV and then check that all 
its elements are actually extreme solutions. 

The coincidence of boundaries is not specific for £, rather it holds for other 
number triangles and more sophisticated graded graphs [13] . A common 
reason for this phenomenon is some law of large numbers, like the law of 
large numbers for exchangeable — 1 random variables in the case of Pascal 
triangle. On the other hand, there are simple examples of graded graphs for 
which the Martin boundary is strictly larger than the extreme boundary [7] . 

3 D— arrangements 

Let [n] = {l,...,n} and let Perm(n) be the set of permutations of [n]. 
We write permutations 7r G Perm(n) in the conventional one-row notation 
7i = 7r(l) . . .7i(n) (and ignore the group structure on Perm(n)). A position 
j G [n— 1] is said to be a descent of n if n(j) > n(j + 1). By D{tx) we denote 
the total number of descents of tt. For instance, tt = 7356241 e Perm(7) has 
descents at positions j = 1,4,6, hence D(7356241) = 3. 

According to a well-known combinatorial interpretation, the Eulerian 

3 The definition corresponds to the entrance boundary in [S]. 




of V. 
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numbers count permutations with a given number of descents: 



Card{7r G Perm(n) : D(n) = k} = Q) , (14) 

as is easily shown by checking that the counts satisfy the recursion (JTJ (or 
see (HI Section 6.2]). We establish next a more delicate connection. 

Observe that removing n from a permutation of [n] yields a projection 
p n : Perm(n) — > Perm(n — 1). For instance, 3412 G Perm(4) is projected to 
312 G Perm(3). Clearly, the preimage of any permutation 7r G Perm(n — 1) 
by p n consists of exactly n elements. 

Lemma 5. There exists a bijection b n between Perm(n) and standard paths 
in the graph £ of length n with the following property: the path b n (7i n ) cor- 
responding to Ti n G Perm(n) passes through the vertices (m, D(jr m )) G £ m 
(m = 1, . . . , n), where is n -\ = p n {^n), ■ ■ ■ , tti = P2{^2) are the iterated projec- 
tions ofir n . 

Proof. Choose 7r n G Perm(ra) and let k = D(-K n ). It is readily checked that 
the preimage Pn+ii^n) C Perm(ra + 1) consists of k + 1 permutations with k 
descents and of n — k permutations with k + 1 descents. Observe that k + 1 
is the number of edges linking (n, k) to {n + 1, k) while n — k is the number 
of edges linking (n, k) to (n + 1, k + 1). It follows that if the desired bijection 
exists for some n then it can be further extended to a bijection for n + 1. 
The assertion follows by induction. □ 



The bijections b n are in no way canonical, because we do not distinguish 
among the edges linking adjacent vertices in £. Still, the way we introduced 
6 n 's takes care of consistency for all n. Indeed, let t n be the operation of 
cutting off the last link in a standard path in £ of length n. Thus t n projects 
standard paths of length n onto standard paths of length n — 1. The consis- 
tency of 6 n 's amounts to the commutation relation t n o b n = b n -\ o p n , which 
holds for all n > 2. 

Let 21 = limPerm(ri) be the inverse limit 4 of the finite permutation 
spaces Perm(n) with respect to p n 's. Elements of 21 are infinite sequences 
(ir n ) of consistent permutations n n G Perm(ra) (n = 1,2, . . .), meaning that, 
for each n > 2, p n (ir n ) = 7i n -i- in extension of LemmaElwe have the following 
corollary. 

4 Another inverse limit, the space of virtual permutations, appears in |14j . 
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Corollary 6. The consistent sequence of bijections (b n ) defines a bijection 
between 21 and the set of infinite standard paths in 8 . The bijection has the 
property that the path corresponding to {ji n ) G 21 passes through the vertices 
(n,D(ir n )), n = 1,2, . . .. 

With each permutation 7r G Perm(n) we associate a total order < on 
the set [n], in which 7r(l) < • • • < ir(n). Likewise, every element (7r n ) G 21 
determines a total order on the set N = {1,2,...} such that, for each n, the 
total order restricted to the subset [n] C N is the one given by 7i n . Conversely, 
any total order on N can be obtained in this way, from some element of 21. 
For this reason, we call the elements of 21 arrangements and identify them 
with the total orders on N.Two obvious examples of arrangements are the 
standard order 1 < 2 < 3 < . . . and the inverse order • • • < 3 < 2 < 1; the 
corresponding paths in £ go along the left side of the Euler triangle and 
along its right side, respectively. 

As a projective limit of finite sets, 21 is a compact topological space. 
Given a probability measure P on 21 we can speak of a random arrangement 
II = (Il n ), where Il n G Perm(n) are consistent random permutations, such 
that the law P n of U n is the pushforward of P by the canonical projection 
21 — > Perm(n). Conversely, by Kolmogorov's measure extension theorem, 
each sequence of distributions (P n ) determines a unique random arrangement, 
provided the sequence is consistent with respect to all projections p n . 

The random arrangements relevant to our discussion have one special 
property of sufficiency. 

Definition 7. We say that a random arrangement II = (Il n ) is a D-arrange- 
ment if for every n — 1, 2 . . . and -n n G Perm(n) the probability of the event 
n n = 7r n depends on the couple (n,D(n n )) only. 

That is to say, for a D-arrangement II = (IT n ), the number of descents is a 
sufficient statistic: conditionally given D(II n ) = k the distribution of II n is 
uniform on the set of permutations of [n] with k descents, for each n and k. 

Two trivial examples of D-arrangements are the nonrandom arrange- 
ments given by the standard order and the inverse order. The corresponding 
measures on 21 are the Dirac masses at points (7r n ) = (l...n) G 21 and 
{^n) = (n . . . 1) G 21, respectively. Notice that these two are the only Dirac 
measures on 21 corresponding to D-arrangements. 
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More substantial example is the random arrangement IT for which ev- 
ery II n has uniform distribution on Perm(n). This is the only exchangeable 
random arrangement, whose probability law is invariant under arbitrary per- 
mutations of the set N. 

Now, Corollary |U] implies: 

Lemma 8. The formula 



defines an affine isomorphism V = (V n k) <-> P = (P n ) between V and the set 
of probability laws for D- arrangements. 

Equivalently, in terms of quantities V n k and random paths in £ corresponding 
to D-arrangements, V n k is the probability that a random infinite path (with 
distribution P) will pass through the vertex (n, k). In the sequel we will not 
distinguish between solutions to (J3J) and random D-arrangements. 

4 Bucket sorting 

Here we use the correspondence of Lemma |H1 for constructing a family of 
solutions V G V. The following algorithm, called bucket sorting, exploits 
a multinomial distribution and is a simplest of the algorithms of this kind, 
widely known in computer science JH] and dynamical systems j2J E3 • 

Fix x G {0,1,2,...} and imagine x + 1 buckets arranged in some order. 
Suppose each of the numbers 1, 2, • • • G N is sent to one of the buckets 
with equal probabilities (x+ l) -1 , . . . , (x+ independently of the other 
numbers. For each n this yields a random allocation of integers 1, . . . , n in 
the buckets. Arranging the integers within each bucket in increasing order 
and putting the resulting sequences together (in the order of the buckets), 
the allocation of n integers is transformed into a random permutation of 
[n]. By the construction, has at most x descents. 

Lemma 9. The infinite sequence = (IL^) produced by the bucket sorting 
is a D- arrangement. The corresponding array = iW^k) G V given by 
formula 



Pn{n n ) = Ki,Z)(7r n )> 



n = 1,2,.. 



7f n G Perm(n) 




(15) 
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Proof. By the very construction, the random permutations are consistent 
with respect to the projections p n , hence is a random arrangement. Given 
7r n G Perm(n), let us compute the probability of the event IT^ = n n . The total 
number of possible allocations of 1, . . . , n into buckets equals (x + l) n , and 
all of them are equally likely. Thus, it suffices to compute the number of the 
allocations resulting in 7r n . Any such allocation is determined by a partition 
of the sequence n n = 7T n (l) . . . 7i n (n) into x + 1 consecutive fragments (some 
of which can be empty), and any such partition can be encoded by placing 
x vertical bars separating the fragments. Observe that for each descent j G 
[n— 1] of the permutation n n at least one bar has to be placed between 7r n (J) 
and 7T n (j + 1). For k = D(jr n ) we see that k positions of bars are fixed by the 
descents, so that the allocation is actually determined by the remaining x— k 
bars, which can be placed arbitrarily. Since the bars are indistinguishable, 
the number of possibilities equals ( n+ ^ _fc )- Thus, the probability of ir n is 
given by the right-hand side of f|15j) . Since this expression depends only 
on k = D(ir n ), we conclude that is a D-arrangement and (fT5|) is the 
corresponding element of V. □ 

Remark 10. The fact that formula ()15j) determines a solution to Q amounts 
to a binomial identity which is easy to check directly: 

„+i )+<»-*)( „ +1 J- 

while the total probability rule ()13|) becomes 

(-^)-=EO( n+ r*)- 

k=0 v 7 

which is equivalent to Worpitzky's identity [8, (6.37)]. 

Recall that in Section El we introduced arrays y N>c = (V^) G 
solving the "iV-truncated" version of recursion (j2J). 

Lemma 11. Fix x G {0, 1, . . . }, let = (11^) be the D-arrangement re- 
sulting from the bucket sorting, and let = {W* k ) G V stand for the 
corresponding array. Then V converge to W* , that is 

lim V£" = W£ (n=l,2,...;0<fc<n-l). 

N—>oo 



11 



Proof. It is more convenient to deal with quantities 

For fixed n, the vectors 

(^•••>KSli) and (W«,..., W- n _ x ) 

are the distributions at time n of the iV-step Markov chain (introduced in 
Section |2J) whose initial distribution at time N is 

(V^,...,V^) or (^ ,..,^.i)- ( 16 ) 

respectively. 

Recall that both V^k and W 7 ^, vanish for k > x and, moreover, 

V& X = <U- (17) 
We claim that it suffices to prove the limit relation 

lim W£ x = 1. (18) 

Indeed, since are probability distributions, (fT7j) and (|T5j) imply that the 
total variance distance between them goes to as iV — > oo, which implies 
the assertion of the lemma. 

To prove (fTHj) we turn to and observe that is just the probability 
for the random permutation LT^ to have the maximal possible number of 
descents x. In terms of the random allocation of the numbers 1, . . . , N, this 
means that all buckets are nonempty and the largest number in each bucket 
(except the last bucket) is larger than the smallest number in the next bucket. 
If this were not the case, all the numbers in one of the buckets were smaller 
than those in the next bucket. Elementary estimates which we postpone to 
the proof of Lemma show that the probability of such an event tends to 
as iV -> oo, which yields (JIBjl. □ 

Remark 12. By virtue of the explicit formula (|15|h the limit relation (fTHj) 
is equivalent to an asymptotic relation for the Eulerian numbers: 

lim ( ) ~ (x+l) N for fixed x = 0, 1, 



N—yoo \ X 



This relation can be readily checked directly. For instance, it follows from 
formula (J22|) below. 
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Lemma ITT1 shows that the family {W*} is contained in the Martin bound- 
ary. Actually, a stronger claim holds: all W*"s are extreme. We show this in 
Lemma El below. But first we will prove the law of large numbers for U*. 

Lemma 13. Fix x = {0, 1, . . . }. We have 

lim D(U*) = x with probability 1. 

n— >oo 

Moreover, this property of is characteristic. 

Proof. Suppose there are just two buckets, x = 2. Then -D(n^) < 2 means 
that for some m G [n] the integers l,...,m fall in the first bucket, and 
m + 1, . . . , n in the second, which is an event of probability (n + l)/2 n . Since 
the series of these probabilities converges, the Borel-Cantelli lemma yields 
the claim. The general case x > 2 is reduced to the estimate in the case 
x = 2 by focussing on two consecutive buckets and using elementary large 
deviation bounds for Bernoulli trials to show that the chance for less than, 
say, nj h integers in both buckets goes to exponentially fast with n. The 
uniqueness follows as in Lemma 1771 □ 

Lemma 14. Elements W* G V resulting from the bucket sorting are extreme. 

Proof. If W* is a mixture of some V 1 , V 2 G V then by the first assertion 
of Lemma 1731 the arrangements corresponding to V\ and V<i must satisfy the 
same law of large numbers as W*. But then by the second assertion of the 
lemma V1 — V2 — W*. Hence is extreme. □ 



5 Proofs of the main results 

We start with reducing the set of parameters needed to determine a generic 
solution V G V. 

Lemma 15. The sequence (V n o) uniquely determines V G V. 
Proof. Writing (j^J as 

1 k + 1 

K+i,fc+i = r V nk V n+hk (0 < k < n - 1), 

n — k n — k 

we see that for each k G {0, 1, ... } the sequence (V^+i : n = k+2, k+3, . . . ) 
is uniquely determined by the sequence {V n k : n = k+1, k+2, . . .). Induction 
in k ends the proof. □ 
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Consider the truncated arrays V N " = {Vjg") e V {N) with parameters N 
and x, < x < iV — 1, introduced in Section 0J As in the proof of Lemma 
El we introduce the modified array V Nlt = (V* k *) with = { n k )V* k *, 
and we recall that the nth row of V N>C is the distribution at time n of the 
Markov chain started at time N from the vertex (N, x). Since V^q = V^*, 
the quantity Vj^* equals the probability of the event that the Markov chain 
will pass through vertex (n, 0). 

Furthermore, there is a monotonicity property analogous to that of gen- 
eralized Stirling triangles in [7j. 

Lemma 16. For fixed n G {1, . . . , iV}, the coordinate V^q* does not increase 
as x varies from to N — 1 . 

Proof. We employ a simple coupling argument. Given two numbers x < x 7 , 
consider two Markov chains which start from vertices (N,*c) and (iV, x'), 
respectively. We settle both chains on a common probability space assuming 
that the jumps are independent as long as the trajectory of the first chain 
does not intersect the trajectory of the second chain, but once the trajectories 
meet, they merge. The merge does not affect the marginal law of each of the 
chains, since both are directed by the same transition probabilities. The key 
property of the coupling is that each trajectory of the first chain remains on 
the left of the trajectory of the second chain, before the trajectories merge. 
Observe now that after reaching the left side of the Euler triangle, a trajectory 
can only process along this side. Consequently, if a trajectory of the second 
chain passes through (n, 0) then the trajectory of the first chain reaches the 
left side of the triangle at some time n' > n, hence passes through (n, 0), too. 
Thus, the chance for the first chain to pass through (n, 0) is not less than that 
for the second chain. This proves the desired inequality > . □ 

We proceed with the proof of Theorems ^ and 121 Our strategy is to 
determine first the Martin boundary by directly identifying all solutions W 
that appear as limits of arbitrary sequences of the form V N '^ N \ 

Assume first that x(N) = x with some fixed x, for all N large enough. 
Then, by Lemma ^2 the sequence V N,y ^- N ' converges to the array given 
by formula (|15|) which we also display here for reader's convenience: 

"S = (" + r*)/ (x+1) " (19) 
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Next, assume that x(iV) = N — 1 — x, where x is fixed. Observe that this 
limit regime is reduced to the preceding one by application of the symmetry 
(n, k) — » (n, n — 1 — k) of the Euler triangle S. Therefore, in this case the 
sequence converges to the array W 1 * with components 

fiS=(* + * + 1 )/<*+ir. (20) 

The random D-arrangement corresponding to can be produced by the 
obvious analogue of bucket sorting in which integers within each bucket are 
arranged in decreasing order. 

Further on, from (|19|) and (J2(Jj) it is readily seen that there exists the limit 

lim W* = lim W" = W°° 

M — >0O X — >00 

with components 

W nk = — y 

Clearly, W°° E V. 

Now we claim that if both x(iV) and N — 1 — x(iV) go to infinity then 
the sequence V N,>t ^ N ' converges to W°°. To that end, observe that a general 
bound 

/ n \ — 1 

o<K fe <( A; ) , (21) 

which follows from Lemma EJ holds for all V G V and implies that V is 
compact in the product topology. By the compactness, passing if necessary 
to a subsequence of (V N '^ N ' > ) we can always achieve convergence to some W, 
hence it is enough to show that W = W°°, and by Lemma ITBl this is further 
reduced to showing that W n0 = = For any fixed x we have 

x < x(N) < N - 1 - x 

for large N. Applying Lemma El we obtain the bound 

> W n0 > W: o (n = l,2,...)- 

Now, sending x to infinity we conclude that W„o = -7, as wanted. 
We have shown that the Martin boundary consists of the elements 

W* (x = 0,1,2,...), (x = 0,1,2,...), and W°°. 
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Comparison with formulas of Theorem ^ shows that these are exactly the 
arrays W{6) with 9 > ~, 9 < |, and 9 = ~, respectively. A remarkable fact 
emerges: the single entry (2, 0) distinguishes all these arrays. 

By Lemma ITU the arrays W 3 * are extreme. By symmetry, the arrays W* 
are extreme, too. To finish the proof of the theorems it remains to check that 
W°° = W{\) is extreme. We postpone this to the next section. 

The Corollary |3] follows from general results. Since the space V is com- 
pact, metrizable and separable, the well-known Choquet theorem [2H1 §3] 
implies that each solution V G V may be represented as a convex mixture of 
the extreme solutions W G ex(V). A simple general argument shows that V 
is a Choquet simplex (that is, the cone generated by V is a lattice cone), see, 
e.g., [T%1 Lemma 9.3]. The uniqueness of representation now follows from 
another Choquet's theorem, see [20J §9]. 

6 End of proof: extremality of the exchange- 
able arrangement 

In this section P denotes the probability measure on 21 corresponding to 
the array W°° = G V. The characteristic property of P is that, for 

each n, the image P n under the natural projection 21 — > Perm(n) is the 
uniform measure assigning to all permutations ir n G Perm(n) equal weights 
-=f. Let II be the random arrangement with law P. This is the exchangeable 
random arrangement, invariant under the natural action on the space 21 of 
permutations of N. 

The following useful construction of II is found in [TJ. Let A 1; X 2 , . . . be 
independent random variables, with uniform distribution on the unit interval 
[0, 1]. The Aj's are pairwise distinct with probability one. Define a random 
total order on N by the rule i < j if < Xj. Clearly, for each n, the 
resulting random permutation Il n of [n] depends only of X\, . . . , X n and, by 
exchangeability of X/s, U n is uniformly distributed on Perm(n). 

Known moments of D(U n ) follow easily from this realization. 

Lemma 17. Let U n be the uniform random permutation of[n}. The random 
variable D(U n ) has mean (n — l)/2 and variance (n — 1)/12. 

Proof. Clearly, D(U n ) equals the number of descents in the random sequence 
Xi, . . . , X n , that is, the number of indices i G [n — 1] such that Xj > X i+ \. 
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Therefore, denoting by Xi the indicator of the event Xi > X i+ % we have 

n-1 
i=l 

The result easily follows from the relations 

E( Xi )=M(Xi) = h E(xiXi+i) = h E( XiXj )=0 (\i-j\>2). 

□ 

Since both the mean and the variance exhibit a linear growth, standard 
application of Chebyshev's inequality gives: 

Corollary 18. Under the uniform distribution, D(U n )/(n— 1) — > | in prob- 
ability. 

We have now all tools to finish the argument of Section |3] by showing 
that W°° is extreme. Assume the contrary, then the boundary ex(V) reduces 
to {W*} U {W 3 *} and hence P can be written as a convex combination of 

the measures P* and P* (the laws of W* and W*), x = 0, 1,2, By 

Corollary there exists a sequence of numbers n\ < < ■ ■ ■ such that 
D(Y[ n .) I \n — 1) — > | almost surely. On the other hand the same ratio goes 

to or 1 under the distribution P* or P 1 *, respectively. This leads to a 
contradiction, so the proof is complete. 

7 Concluding remarks 

7.1 Permutations with descent— set statistic 

We were led to consider the Eulerian triangle £ in connection with a larger 
graded graph Z of zigzag diagrams 0. With edge multiplicities taken into 
account, both graphs have the same path spaces, but Z has more vertices 
and much more rich branching. The boundary problem for Z amounts to 
describing all random arrangements II = (II n ) with the property that the 
distribution of each Il n is uniform conditionally given the set of descent posi- 
tions in II n . In jU] we established that the distribution of a random total order 
determined by such II must be spreadable, that is invariant under increasing 
mappings N — > N. D-arrangements are the simplest of this kind, and the 
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extreme D-arrangements we described here are also extreme solutions to the 
boundary problem for Z. 

Analogous connection exists between Kingman's graph K, of partitions 
and the Stirling triangle S of the first kind [7 . The relevant random objects 
are exchangeable partitions of N and a smaller class of partitions which have 
the number of blocks as sufficient statistic. In that case the situation is more 
interesting than the one for Z and £: extremes solutions to the boundary 
problem for S (the celebrated Ewens partition structures) are decomposable 
along the boundary of /C, with the mixing measure being the remarkable 
Poisson-Dirichlet distribution [2*2*] . 

7.2 A problem of moments 

Corollary El and (JJJ) tell us that every solution V G V satisfies 



for some unique probability distribution p on the parameter set G. An inverse 
problem asks one to characterize all sequences (V n o) with Viq = 1 which can 
be represented in this form. An answer is suggested by the argument in 
Lemma [T3] which says that there is a linear operator V : (V n o) i— > (V n k) 
which maps an arbitrary sequence to a solution of So the necessary and 
sufficient condition for representability is that V applied to (V^o) produces a 
nonnegative array. 

The analogous question for Pascal's triangle is the Hausdorff moment 
problem, with kernel 8 n , where 9 ranges in [0,1]. In this classical case the 
analogue of V associates with each sequence the array of its iterated differ- 
ences, whose positivity is Hausdorff's condition called total monotonicity. 

7.3 Remarks on the uniform case 

The Eulerian numbers are given by the formula 





(22) 
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which compared with Laplace's formula in [5J Section 1.9] shows that 

Prob(A; < Y x + . . . + Y n < k + 1) = Q) /n! 

for Yi, Y2, . . . independent random variables with uniform distribution on 
[0, 1]. The following explanation of this coincidence is borrowed from p. 
296]. For x > let x = [x\ + {x} be the decomposition of x into integer and 
fractional parts, with < {x} < 1. Consider 

S J = Y 1 + :. + Y j , X, = {5,}, 

then Xi, X 2 , . . . are also independent, uniform on [0, 1]. Observe that Xj > 
Xj + i each time [Sj\ < [Sj + i\ = \_Sj\ + 1 and recall the discussion preceding 
Lemma El 

Improving upon Corollary we see that the convergence D(U n )/(n — 
1) — > I holds with probability 1. The central limit theorem applied to 5„'s 
entails that the distribution of D(U n ) is asymptotically Gaussian. This con- 
nection with sums of random variables has been a starting point for many fine 
results on descents in uniform permutation. See for recent development 
and references. 
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